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Based on the analysis of data from the numerous dedicated experiments on plasma dis¬ 
ruptions in the TEXTOR tokamak the mechanisms of the formation of runaway elec¬ 
tron beams and their losses are proposed. The plasma disruption is caused by strong 
stochastic magnetic field formed due to nonlinearly excited low-mode number magneto- 
hydro-dynamics (MHD) modes. It is hypothesized that the runaway electron beam is 
formed in the central plasma region confined by an intact magnetic surface due to the 
acceleration of electrons by the inductive toroidal electric field. In the case of plasmas 
with the safety factor q{ 0) < 1 the most stable runaway electron beams are formed by the 
intact magnetic surface located between the magnetic surface q = 1 and the closest low- 
order rational surface q = m/n > 1 ( q = 5/4, q = 4/3, ...). The thermal quench time 
caused by the fast electron transport in a stochastic magnetic field is calculated using 
the collisional transport model. The current quench stage is due to the particle transport 
in a stochastic magnetic field. The runaway electron beam current is modeled as a sum 
of toroidally symmetric part and a small amplitude helical current with a predominant 
m/n = 1/1 component. The runaway electrons are lost due to two effects: (i) by outward 
drift of electrons in a toroidal electric field until they touch wall and {ii) by the forma¬ 
tion of stochastic layer of runaway electrons at the beam edge. Such a stochastic layer 
for high-energy runaway electrons is formed in the presence of the m/n = 1/1 MHD 
mode. It has a mixed topological structure with a stochastic region open to wall. The 
effect of external resonant magnetic perturbations on runaway electron loss is discussed. 
A possible cause of the sudden MHD signals accompanied by runaway electron bursts 
is explained by the redistribution of runaway current during the resonant interaction of 
high-energetic electron orbits with the m/n = 1/1 MHD mode. 


1. Introduction 

One of the severe consequences of the plasma disruptions in tokamaks is the generations 


of the runaway electron (RE) beams (see, e.g., Wesson et al. (19891; Gill (1993); Schuller 
fll995l ); |Gill et al.\ pOOP) [2002] ); |Wesso^] fl2004D ; |Boozer| fl2012| ); |Papp et al\ ( |2013D and 
references therein). The REs generated during the disruptions of tokamak plasmas may 
reach several tens of MeV and may contribute to the significant part of post-disruption 
plasma current. The prevention of such RE beams is of a paramount importance in future 
tokamaks, especially in the ITER operation, since it may severely damage a device wall 
(Beco ulet et q(.| [2013). 
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At present there are several proposals to mitigate REs generated during plasma dis¬ 
ruptions. The mitigation of REs by the gas injections has been discussed (see, e.g., Refs. 


(Hender et al. 2007 Whyte et al. 2002 2003 Bakhtiari et al. 2002 2005 Granetz et al. 

2007 

Bozhenkov et al. 2008 Pautasso et al. 2009 

Lehnen et al. 

2009 

Hollmann et al. 

|2010 

Reux et al: 2010j Lehnen et al. 20111). Suppres 

sion of REs by the resonant magnetic 
ussed since late 1990s (see, e.g., Refs. 

perturbations (RMPs) has been also intensively disc 

(Kawano et al. 1997[ Tokuda & Yoshino 1999 Helander et al. 2000 Yoshino & Tokuda 


up to now there is no a regular strategy to solve this problem. One of the reasons is 
that the physical mechanisms of the formation of REs during plasma disruptions is still 
not well- known. The different scenarios of runaway formation during plasma disruptions 


are discussed in literature. Particularly, in Refs. (Frilop et al. 2009 Fiilop & Newton 


2014) the possible roles of whistler waves on the generation of REs and Alfvenic wave 


instabilities driven by REs have been discussed. 

There were numerous dedicated experiments to study the problem of runaway current 


generation during plasma disruptions triggered by massive gas injections (MGI) in the 

TEXTOR tokamak (see, e.g., ( 

Forster et al. 

2012 Zeng et al. 

2013 

Wongrach et al. 

2014)), in KSTAR tokamak (Chen et al. 2013 

, the JET tokamak 

(Plyusnin et al. 2006 

Lehnen et aZ.||2011), in DIII-D 

Hollmann et al. 

2010 Commaux et al. 

2011 Hollmann 

et al. 2013), Alcator C-Mod (Olynyk et al. 2013 

, and others. In these works the depen- 


dencies of RE generation on the toroidal magnetic field, on the magnetic field fluctuations, 
on the species of injection gases have been investigated. Particularly, in KSTAR tokamak 
(Chen et q/.|2013) it has been found that there is no the toroidal magnetic field threshold 


Bt < 2 T as was indicated by previous experiments in other tokamaks. In Ref. (Izzo 
et aZ.|[20lI 2012) MHD simulations have been performed to study the confinement REs 


generated during rapid disruptions by MCI in DIII-D, Alcator C-Mod, and ITER. Such 


simulations with two different MHD codes have been carried out by Izzo et al. (2012) to 


analyze shot-to-shot variability of RE currents in DIII-D tokamak discharges. 

These numerous experiments show the complex nature of plasma disruption processes 
especially the formation of RE beams, and its evolution. One of the important features 
of this event is its irregularity and variability of RE beam parameters from one discharge 
to another one. This indicates the sensitivity of disruption processes and RE beam for¬ 
mations on initial conditions which is the characteristic feature of nonlinear processes, 
particularly, the deterministic chaotic system. Therefore, ab initio numerical simulations 
of these processes may be not always successful to understand their mechanisms because 


of the complexity of computer simulations of nonlinear processes (Kadanoff 2004). The 


problems of numerical simulations of plasma disruptions is comprehensively discussed by 


Boozer (2012). 


In this work we intend to approach to this problem from the point of view of Hamil¬ 
tonian chaotic systems, mainly the magnetic stochasticity in a magnetically confined 


plasmas (Abdullaev 2014). Based on the ideas of these systems and analyses of numer¬ 


ous experimental results, mainly obtained in the TEXTOR tokamak we propose possible 
mechanisms of formation and evolution of RE beams created during plasma disruptions. 
Since a self-consistent theoretical treatment of all these processes is very complicated we 
developed theoretical models for each stages of a plasma disruption. These models are 
used to estimate the characteristic times of the thermal and current quenches, the spatial 
size of runaway plasma beam and their decay times, the speed of RE radial drifts, the 
effect of magnetic perturbations. 

It is believed that the plasma disruption starts due to a large-scale magnetic stochas¬ 
ticity caused by excited of MHD modes with low poloidal m and toroidal n numbers, 
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(m/n = 1/1,2/1,3/2, 5/2,...) (see, e.g., (Wesson 2004 
[Schuller 1995J and references therein). The heat and particle transports in the strongly 
chaotic magnetic field causes the fast temperature drop and ceases the plasma current. 
However, at the certain spectrum of magnetic perturbations, for example, at the suffi¬ 
ciently small amplitude of the m/n = 1/1 mode the chaotic field lines may not extend to 
the central plasma region due to the creation of an intact magnetic surface. In the case 
of plasmas with the safety factor q{ 0) < 1 at magnetic axis p = 0 the intact magnetic 
surface can be located between magnetic surface q = 1 and the nearest low-order rational 
surface q = 5/4 [or q = 4/3, ...]. This intact magnetic surface confines particles in the 
central plasma region and serves as a transport barrier to particles during the current 
quench. Electrons in the confined region are accelerated due to large toroidal electric 
field and forms the relatively stable of RE beams. 

This occurs, for instance, when the plasma disruption initiated by the heavy Ar gas 
injection which does not penetrate deep into the plasma, therefore it does not excite the 
m/n = 1/1 mode with the sufficiently large amplitude. In contrary, the injection of the 
lighter noble gases neon and helium does not generate runaways. The reason is that light 
gases penetrate deeper into the plasma and excite the large-amplitude (m/n = 1/1) 
mode. 

The existence of an intact magnetic surface and its location depends on the radial pro¬ 
file of the safety factor and the spectrum of magnetic perturbations. The latter sensitively 
depend on the plasma disruption conditions and vary unpredictably from one discharge 
to another during plasma disruptions. This makes RE formation process unpredictable 
and may explain a shot-to-shot variability of the parameters of RE beams. 

The role of the safety factor profile in the formation of RE beams can be pronounced 
during disruptions of plasmas with the reversed magnetic shear. In the plasmas with 
the non-monotonic radial profiles of the safety factor it has been observed an improved 
confinement of energy and particles due to the internal transport barrier located near 


the minimal value of the safety factor, i.e., near the shearless magnetic surface (Levinton 
et aL|199~5 Strait et ~o77|fl995 ). During disruptions this magnetic surface acts as a robust 


magnetic barrier that separates a chaotic magnetic field formed in outer region from the 
penetration into the central plasma region. Electrons confined by the shearless magnetic 
surface can form a stable RE beam with a relatively large transversal size. Recently 
published results of the disruption experiments in the TFTR tokanrak with the reversed 
magnetic shear indeed show the formation of a large RE beam with long confinement 
times ( Fredrickson et aZ.||2015 >. 

Based on this mechanism we study the main three stages of the post-disruption plasma 
evolution: the fast thermal quench (TQ), the current quench (CQ), and the RE beam 
evolution. The physical processes during each of these stages will be studied by theo¬ 
retical models. These processes are the formation of stochastic magnetic field, heat and 
particle transport in a stochastic magnetic field, the acceleration of electrons by induc¬ 
tive electric field, the lost mechanisms of REs, and the effect of internal and external 
magnetic perturbations. The short report on this study is to publish in Abdullaev et al. 
p015j ). 

The paper consists of eight sections. Mathematical tools and models employed to study 
the problems are given in Supplementary part. The numerous data obtained during 
the dedicated experiments in the TEXTOR tokamak are analyzed in Sec. [2] Possible 
mechanisms of plasma disruptions with a RE beam formation is proposed and analyzed in 
Sec. [3} The transport of heat and particles during the fast TQ and the CQ stages of plasma 
disruption are studied in Sec. [4} The model of a post-disruption plasma beam is proposed 
in Sec. [5j Using this model a time-evolution of guiding-center (GC) orbits of electrons 
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accelerating by the inductive toroidal electric field is studied in Sec. [6] Particularly, the 
change of RE confinement conditions with decreasing the plasma current and increasing 
the electron energy and the outward drift of GC orbits are investigated. The effect of 
external and internal magnetic perturbations on the RE confinement are discussed in 
Sec. [7J In the final Sec. [8] we give the summary of obtained results and discuss their 
consequences. 


2. Description of plasma disruptions 

The TEXTOR was a middle size limiter tokamak with the major radius Rq = 1.75 
m, the minor radius a = 0.46 cm. The toroidal field Bq can be varied up to 2.8 T, and 
the plasma current take up to 600 kA. In the experiments the plasma disruptions were 
triggered in a controlled way by gas injections using a fast disruption mitigation valve 


(DMV) ( Bozhenkov et «/.|2007 Finken et al. 2008| 2011| Bozhenkov et at|2011 ). Particu¬ 
larly, the disruptions with REs were triggered by argon (Ar) injection. The runaway-free 
disruptions were triggered either by helium or neon (He/Ne) injection performed by the 
smaller valve. The effect of the externally applied RMPs on the REs generations has 
been investigated using the dynamic ergodic divertor (DED) installed in the TEXTOR 
tokamak. 

Below we analyze the experimental results of discharges with the predisruption plasma 
current I p = 350 kA and the toroidal field B t = 2.4 T. Figure [l] (a) illustrates typical 
disruptions of the discharges of the TEXTOR tokamak with and without RE generations. 
Specifically, it shows the time evolution of plasma parameters (the loop voltage Vi oop , 
the electron cyclotron emission (ECE), the soft X-ray (SXR) signal, the Mirnov signal, 
and scintillation probe signals (ScProbe) during disruptions of the discharges with REs 
(#117434, #117859, #119978, #120140) and without (#117444) REs. 

There are also some discharges with untypical RE currents and shorter current decay 
times. The two examples of such discharges are shown in Fig.[l](b). We will discuss some 
features of these discharges at the end of the section. 

The typical behavior of the plasma during the disruptions is following. The gas (Ar or 
Ne/He) was injected at the time instant t = 2 s. One can distinguish three stages of the 
disruption with the REs: the first (or fast) stage in which a sudden temperature drop 
occurs, in the second stage the plasma current starts to decay with a higher rate, and 
in the third stage the current decay slows down and the current beam with the REs is 
formed. 


2.1. Thermal quench stage 

The first fast stage starts after a few milliseconds (between 2 ms and 5 ms) after the gas 
injection and ends with a sudden temperature drop (a thermal quench) in a time interval 
about one ms as seen from the ECE signals shown in a detail in Fig. [l] (c). The Mirnov 
signals indicating magnetic activities start just before of this time interval and they last a 
few milliseconds until a significant decay of the plasma current for the RE-free discharges 
or establishing the current with the REs (see Fig. [I] (a)). The close-up views of the ECE 
signals and the Mirnov signals in this stage for the discharges with REs (#117434) and 
without REs (#117444) are shown in Figs. [l] (c) and (d), respectively. 

For our study it is of importance to analyze in details the difference in the time 
development of the TQ stage of disruptions without REs and with REs, initiated by the 
massive injection of lighter (Ne) and heavier (Ar) noble gases, (i) In Fig. [l] (c) one can 
see that the TQ stars roughly 2.6 fis, in the former case, and 3.7 fi s, in the latter one, 
after the initiation of injection. The ratio of these delay times is of 0.7 « l/y/2 and can 
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Figure 1. (a) Time evolution of the disruption of the TEXTOR shots #119978 (black solid 
lines) and #117444 (red curves) (from top to bottom): the plasma current, the loop voltage, 
the ECE signal, the SXR signal, the Mirnov signal, and scintillation probe (Sc.Probe) signal, 
(b) The same but for the discharges #117859 (blue curves) and #120140 (black curves), (c) 
Initial stage of the temporal evolution of the plasma current (solid curve on the l.h.s. axis), 
ECE signals (r.h.s. axis); (d) the Mirnov signals (r.h.s. axis) during a plasma disruption with 
(#117434 and #117507) and without (#117444) RE generations. Ip RE ^ is the initial value of 
the plasma current with REs. Disruptions for discharges #117434, #119978, #117507, #117859 
#120140 are initiated by Ar injections, and #117444 by Ne injections. 
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Figure 2. Dependence of the penetration length l g on the electron temperature T e for He, Ne 
and Ar atoms. The electron density n e = 10 19 m~ 3 . 


be well explained by the difference in the atomic weight A g of gases in question. Indeed, 
the flow velocity V g ~ 1 / s/A^ of neon atoms is by a factor of \[2 higher than that for 
argon and, thus, atoms enter the plasma after a respectively shorter time. 

(ii) By comparing Figs, [l] (c) and[l](d) we find that by injection of neon the TQ stage 
is finished before MHD perturbations are triggered. That is in this case TQ is completely 
due to cooling induced by the presence of impurity atoms squeezed in a narrow jet and 
penetrating deep enough into the plasma core. The penetration depth of gas atoms, 
l g = V g / ( ki on n e ), where fc 9 on is the ionization rate coefficient. Figure [ 2 ] displays l g versus 
the electron temperature T e with the electron n e = 10 19 m -3 computed for He, Ne and 
Ar, by using open atomic data base ADASjf] for fc 9 on and assuming that the gas jet has 
a radial velocity of 2 sound speeds at the room temperature. One can see that for light 
gases, He and Ne, l g can exceed the minor radius of TEXTOR of 0.46 m if inside the 
gas jet the plasma is cooled down to a temperature of several eVs by energy losses on 
excitation and ionization of gas atoms and thermalization of generated electrons, as it is 
demonstrated in Koltunov & Tokar (20111. The rest of magnetic surfaces is cooled down 


by the heat conduction along magnetic field lines to the jet area (Tokar & Koltunov 


2013). Thus, the plasma is cooled down as a whole during a time of a/V g « 0.5 — 1 ms. 


Only later tearing modes are exited on numerous resonant magnetic surfaces due to the 
growth, as l/T^ 2 , of the plasma resistivity. 

(iii) In the case of Ar the penetration depth l g is much smaller, mostly due to much 
larger ionization rate, than for He and Ne atoms and Ar gas is ionized at the plasma 
edge. Due to edge cooling very sharp radial gradients of the plasma resistivity and plasma 
current density j are generated. Since the growth rate of tearing modes ~ ( dj/dp ) 4 ^ 5 
(Wesson 2004), MHD modes resonant mostly on outer magnetic surfaces with the safety 
factor noticeably large 1 are triggered. The magnetic field stochastization due to these 
modes leads to the fast cooling of the main plasma volume during a time of 0.1 — 0.2 ms 
(see Sec. 4.1). This explains why in shot 117434 with Ar injection the temperature drop 
happens although later but faster than in shot 117444 with Ne, and MHD activity starts 
to develop even before the TQ, see Figs, [l] (c) and (d). 


2.2. Current quench stage 

The second stage of the plasma disruption begins with the current decay within a mil¬ 
lisecond after the TQ. Particularly, for the discharges #117434 and #117444 the current 

f OPEN-ADAS: Atomic Data and Analysis Structure, http://open.adas.ac.uk/. 
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Figure 3. (a) Determination of the average current decay rate (dl p /dt) at the second and the 
third stages of the plasma disruptions. Solid curve describes the measured time-evolution of the 
current I P {t), dashed straight line f(t) = a — bt approximates the average current decay. The 
coefficient b gives the estimation of \(dl p /dt)\. Curve 1 corresponds to the current decay stage, 
and curve 2 corresponds to the RE current decay stage. 


(a) (b) 




I{, RE) [kA] I ( p RE) IkA] 

Figure 4. (a) Dependence of the average current decay rate \ (dl p /dt)\ on the initial RE current 
I P RE ' > in the current quench regime for a several discharges; (b) The same as in (a) but in the 
RE current decay stage. 


decay starts in 0.47 x 10~ 3 s and 0.87 x 10 -3 s, respectively, after the temperature drop 
(see Figs. [I] (a) and (c)). In discharges without the RE formation the current decays with 
the same rate until it completely disappears in a few millisecond. In the discharges with 
the RE formation the strong current decay stops at the certain value of I„ = Ip RE ' > and 
replaced by slower decay. The initial RE current Ip RE ^ is shown in Figs. Ill (a) and (b). 
In this stage the loop voltage starts to rise due to inductive electric field opposing to the 
current decay. 

The time dependence of the plasma current I p in this stage for all discharges is well 
approximated by the linear function of time I p = I p o — bt, where the coefficient b = 
— ( dl p /dt ) determines the average current decay rate. The scheme of determination of 
b = \(dl p /dt)\ is shown in Fig. [3] The values of the current decay rate ( dlp/dt ) during 
the CQ and RE plateau regimes, the initial RE current I p RE ' for a number of discharges 
are listed Table [1] It also shows the time t max when the applied RMPs, i.e., the DED 
current Ided reaches its maximal value, and the toroidal mode n of the RMPs. 

The dependencies of \(dl p /dt)\ on the initial RE current I p RE ^ for a number of dis- 
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Discharge No. 

2— nd stage 

3— rd stage 

I ( P RE) [kA] 

DED tmaxi Ided 

117434 

40.2 ± 0.45 

4.63 ±0.09 

220.1 

no 

117444 

47.8 ±0.21 

N/A 

N/A 

no 

117507 

42.7 ±0.71 

4.29 ± 0.14 

193.9 

no 

117527 

39.4 ±0.29 

5.65 ±0.12 

220.9 

no 

117543 

50.2 ±0.19 

NA 

N/A 

no 

117859 

33.2 ± 1.11 

11.70 ±0.50 

296.0 

no 

119868 

40.05 ±0.31 

2.45 ±0.01 

224.9 

no 

119869 

39.96 ±0.67 

2.20 ± 0.02 

230.0 

2.02 s, n = 1, 1 kA 

119870 

36.5 ±0.96 

2.23 ± 0.01 

254.8 

2.02 s, n = 1, 1.5 kA 

119874 

37.9 ±0.60 

2.25 ±0.01 

243.8 

no 

119877 

45.3 ±0.47 

25.13 ±0.98 

166.8 

1.9 s, n = 1, 2 kA 

119978 

38.8 ± 0.76 

3.38 ±0.02 

218.4 

no 

119990 

41.98 ±0.23 

3.35 ±0.03 

204.0 

no 

120106 

42.8 ±0.42 

2.03 ±0.03 

196.0 

2.0 s, n = 2, 4 kA 

120107 

40.73 ± 0.29 

2.52 ±0.03 

206.0 

2.0 s, n = 2, 4 kA 

120108 

42.71 ±0.32 

2.25 ±0.02 

191.0 

1.9 s, n = 2, 4 kA 

120109 

42.91 ±0.35 

1.93 ±0.02 

176.0 

1.9 s, n = 2, 4 kA 

120123 

36.4 ±0.47 

3.66 ±0.01 

236.0 

no 

120126 

43.3 ± 0.25 

1.92 ±0.03 

176.0 

2.0 s, n = 2, 7 kA 

120134 

45.0 ±0.97 

2.99 ±0.17 

176.0 

2.0 s, n = 2, 7 kA 

120135 

44.3 ± 0.43 

2.20 ± 0.05 

177.0 

2.0 s, n = 2, 7 kA 

120140 

47.5 ± 0.86 

13.91 ±0.35 

165.0 

1.97 s, n = 2, 6 kA 

120141 

42.2 ±0.30 

2.20 ± 0.03 

194.8 

2.07 s, n = 2, 6 kA 


Table 1 . Parameters of discharges: 1-st column — the discharge number; 2-nd and 3-rd columns 
— the average decay rates \(dlp/dt)\ [in MA/s] of the plasma current I v (t ) in the second and 

(RE') 

the third stages; 4-th column — the initial current of the RE beam Ip ; 5-th column shows 
the parameters of the RMPs, a time tmax when the DED current reaches its maximum values 
Ided, the toroidal mode n. Note, that the discharges #417444 and #117543 are RE-free. 


charges are plotted in Fig. [4] (a) and (b) in the CQ regime and the RE current decay 
stage, respectively. The current decay rate | (dl p /dt ) | for all discharges are the same order 
and lies between 32 and 50 [MA/s] as listed in the 2-nd column of Table [l] and shown 
in Fig. [4] (a). The highest values of \(dl p /dt)\ is observed for the discharges without REs 
and with the lowest values of the RE current I p RE \ As can see from Fig. [ 4 ] there is a 
clear regular dependence of the current decay rate on the initial RE current, which can 
be fitted by a linear function, \(dl p /dt)\ ss A — BI p RE ^ with constant parameters A,B. 
In Sec. 4.2 we will discuss the possible mechanism of such a current decay related with 
the transport of particles in a stochastic magnetic field. 


2.3. RE plateau stage 

In the third stage (RE plateau) of the disruption the rapid current decay is replaced by 
it’s slow decay and it starts the formation of the REs due to the acceleration of electrons 
in the inductive toroidal electric field and the secondary generation of REs. The values of 
the current decay rate | ( dl p / dt) | along with the initial values of the plasma current I p RE ' > 
in this stage for several discharges are listed in the 3-rd and the 4-th columns of Table |T] 
and plotted in Fig. |4] (b). The average values of \{dl p /dt)\ for almost all discharges are 
confined in the interval (2.2, 5.6) MA/s, i.e., in one order lower than the current decay 
rate in the second stage. The values of I p RE ^ are also confined in the range between 170 
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kA and 260 kA, except of some untypical discharges, which will be discussed in the next 
subsection. These values of \(dl p /dt)\ and Ip RE ' > are close to the ones observed in the 
similar experiments in the DIII-D tokamak (see, e.g., (H ollmann et al.pO lO 2013)). 

One should also note that in the RE plateau stage at certain time instants one observes 
also a sudden current drop accompanied by magnetic activity and RE bursts as seen from 
Figs.[l](a) and (b) (see also e.g., Refs. ( Gill et aL|2000 Forster et al.|2012 )). These events 
are probably related to the nonlinear interaction of high-energetic electrons with MHD 
modes which leads formation a stochastic layer at the beam edge open to the wall. We 
will discuss this phenomenon in Sec. |7.3| In the final termination stage one observes the 
quick RE current losses accompanied by magnetic activity. 


2.4. Untypical discharges with REs 

As was mentioned above there are several untypical discharges for which the rates 
(dlp/dt) take highest or lowest values (see Figs. [Tl (b), Fig. [4] (b), and Table [l]). Par¬ 
ticularly, the current decay rate (in the 2-nd stage) for #117859 is lowest and highest 
for discharges #119877, #120140. The RE current decay rate (in the 3-rd stage) for 
these discharges takes highest values. The quantity Ip 11 ^ takes the lowest value for the 
discharges #119877, #120140 and the highest value for #117859 as shown in Fig. |4](b). 
One can notice strong spikes in the SXR signals of these discharges in compared to typ¬ 
ical discharges (see Fig. [I] (a) and (b)). Moreover, the above mentioned bursts of REs 
accompanied by magnetic activities are more pronounced in these discharges. We will 
discuss the peculiarity of these discharges in Secs. [3] and [5j 


2.5. Effect of the RMPs on RE generation 

In a number of discharges the effect of the DED of the TEXTOR (see Sec. 5.2 of Supple¬ 
mentary Part) on the RE generation has been studied. It was found that the RMPs do 
not completely eliminate the RE formation, but it can increase the decay rate ( dl p /dt) 
and decrease Ip RE K This effect depends on the operational mode n, the amplitude of 
the DED current Ided , on the time t max when the maximal DED current is reached. 
As seen from Table [I] the maximal effect is obtained when the maximal Ided. is reached 
before the gas injection at t = 2.0 s, i.e., t max < t = 2.0 s. However, at t max > t = 2.0 s 
the RMPs does not affect at al or it is a very weak. Other experimental observations in 
the TEXTOR also confirm these observations (Koslowski et al. 20141. We will discuss 
this problem in Sec. |7.2| 


3. Formation of a confined plasma beam 

3.1. Main conjecture 

It is believed that the plasma disruption is caused by a large scale magnetic stochasticity 


of field lines due to interactions of nonlinearly destabilized MHD modes (Carreras et al. 
Kadomtsev 1984 Lichtenberg||1984| |Fukuyama et «i.||l998[ |Wesson 2004| |Kruger 


1980 


et al. 


2005 


White 2014)). The global stochasticity are mainly due to the interactions 


of coupled MHD modes with low (m,n) numbers: (m = 1,2,...), (n = 1,2...). The 
structure of a stochastic magnetic field mainly depends of the amplitudes B mn of MHD 
modes and the radial profile of the safety factor q(p), where p is the minor radius of a 
magnetic surface. Depending on these parameters the stochastic magnetic field may fill 
entirely plasma region so that the plasma particles are transported out along chaotic 
magnetic field lines which leads to cease of a plasma current. However, at certain condi¬ 
tions the stochastic magnetic field may not extend up to the central plasma region due 
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to the formation of a magnetic barrier by the outermost intact magnetic surface at p c . 
The electrons confined by this magnetic surface are accelerated by the toroidal electric 
field induced the current decay from outer plasma region, and thus forming a RE beam. 
Let I p (p) be the plasma current flowing inside the magnetic surface of radius p, 

I p (p) = 2t t f j(p)pdp, (3.1) 

J 0 

where j ( p ) is the current density. Then the initial RE current Ip RE ' 1 is mainly determined 
by the predisruption plasma current distribution I p (p) confined by the intact magnetic 
surface p c , i.e., I p RE) « / P (p c ). 

As will be shown in Secs. [6] and [7] the decay of the RE beam mainly depends on the two 
effects: the outward drift of RE orbits induced by the toroidal electric field E v and their 
resonant interactions with helical magnetic perturbations. The outward drift velocity Vdr 
is determined of E v and the RE current, 

v dr oc Ey/lW oc E^/pl- (3.2) 

The most stable RE beams are expected to form when the corresponding drift velocity 
is lowest and the low-order rational magnetic surfaces within the RE beam are absent 
or only one. 


3.2. Possible generic structures of stochastic magnetic fields 
Below we study possible structures of the stochastic magnetic field which may lead to the 
formation of the RE beams. We consider the two type of the safety factor profiles of q(p): 
( i ) the monotonic radial profile and (ii) the non-monotonic radial profile, corresponding 
to the plasmas with the reversed magnetic shear. 

The models for the radial profiles of the plasma current J p (p), the safety factor q(p) of 
the pre-disruption equilibrium plasma, and the MHD magnetic perturbations are given 
in Sec. 3 of Supplementary part. The perturbation magnetic field simulating low-mode 
number MHD modes is given by the toroidal component of the vector potential 

A^(R,Z,tp,t) = - —9-^2 m~ 1 a mn (p) cos (mid -mp + f2 m „i), 

ran 

Q'mn(p) = dd mn U mn (p) (3-3) 

with the mode amplitudes B mn and rotation frequencies Here Bq is toroidal field 
strength, Rq is the major radius Rq , and the functions U mn (p) describes the radial profiles 
of modes. 

One should note that the structure of magnetic field lines in the presence of magnetic 
perturbations is less sensitive to the radial profiles of U mn (p). It is mainly determined 
by the safety factor profiles and the mode amplitudes a mn (p) at the resonant surfaces 
p = Pmm q(Pmn) = m/n (see Sec. 7.2 of Supplementary part). 

Monotonic radial profile of q(p): The case q/O) < 1. The typical TEXTOR 
plasma has the monotonic safety factor profile with the value g(0) < 1 at the magnetic 
axis p = 0. In this plasma, the m/n = 1/1 mode should play an important role on the 
structure of stochastic field lines near the plasma center. At low amplitudes of this mode 
the global stochastic field lines may not reach q = 1 magnetic surface and may form a 
confined region about the plasma center where REs can be generated. At high amplitudes 
of the m/n = 1/1 mode the stochastic field lines may cover entire plasma region with no 
confined particles. 

As was mentioned above in TEXTOR experiments plasma disruptions with REs were 
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Figure 5. Poincare sections of magnetic field lines in a pre-disruption plasma caused by several 
MHD modes: (a) the mode amplitudes b mn = Bmn/B 0 are ( 1 , 1 , 1,1 )xe MH D', (b) b mn = Bm^/Bo 
are (0.5,1,1,1 )cmhd- Red curve corresponds to the last intact magnetic surface, blue curves 
are the resonant magnetic surfaces q = 1, q = 4/3, q = 3/2, and q — 2, respectively. The 
dimensionless perturbation parameter cmhd = 1.5 x 10 -4 . The plasma current I p = 0.35 MA, 
the toroidal field Bo = 2.5 T, the safety factor at the magnetic axis is q(0) = 0.8. 


deliberately caused by the injection of Ar gas while the RE-free disruptions are triggered 
by He/Ne injection. Experiments show that the penetration lengths of atoms depends 
on their atomic weights (Bozhenkov et al. 2008): He (or Ne) atoms penetrate deeper 
into plasma than Argon atoms. The injection of these gases may finally give rise to 
different spectra of amplitudes of MHD modes. One can expect that the amplitude of 
the m/n = 1/1 MHD mode excited by the He/Ne injection is higher than in the case of 
Argon gas injection. 

The two possible distinct generic structures of a stochastic magnetic field before the 
current quench with the RE-free discharge and with the RE discharge are shown in 
Figs. [ 5 ] (a) and (b) by the Poincare sections of magnetic field lines. It is assumed that the 
perturbation field contains several MHD modes: ( m/n = 1 / 1 ), {m/n = 2 / 1 ), {m/n = 
3 / 2 ), and {m/n = 5 / 2 ). In the case shown in Fig. [ 5 ] (a) the normalized mode amplitudes 
bmn = B mn /B 0 are ( 1 , 1 , 1 , 1 ) x c M hd, and in Fig. [ 5 ] (b): ( 1 / 2 , 1 , 1 , 1 ) x c M hd- The 
toroidal field magnitude is Bq = 2.5 T and the dimensionless perturbation parameter 
£mhd = 10 ~ 4 . As seen from Fig. [ 5 ] (a) for the larger amplitude of the {m/n = 1 / 1 ) 
mode the stochastic magnetic field extends up to the central plasma region destroying 
the separatrix of the m = n = 1 island. For the low-amplitude of the {m/n = 1 / 1 ) mode 
shown in Fig. [ 5 ] (b) the stochastic magnetic field does not reach the q = 1 magnetic 
surface and covers the region outer the q = 1 magnetic surface. The last intact magnetic 
surface p c (red curve) is located between the resonant surfaces q{pi) = 1 and q{p$) = 4/3 
(blue curves). 

As seen from Fig. [5] particles in the plasma core are confined by intact magnetic 
surfaces located between resonant surfaces q = 4/3 and q = 1. Plasma beam confined in 
this area is relatively stable. It contains only the m/n = 1/1 MHD mode which does not 
lead to a global stochasticity. The radial transport of particles from the confined area 
can take place only due to small-scale turbulent fluctuations and therefore it has much 
smaller rate than those in the stochastic zone. The confinement time of these electrons 
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(b) 




Figure 6. The same as in Fig. [5]but for the case g(0) > 1: (a) q(0) = 1.1; (b) q{ 0) = 1.2. It is 
assumed the magnetic perturbation contains the m/n modes (n = 1 — 3, m = 1 — 8) with the 
equal amplitudes. The dimensionless perturbation parameter 6mhd = 1-0 x 1CF 4 . The plasma 
parameters are the same as in Fig. [5] 


is sufficiently long enough to be accelerated by the inductive electric field, thus creating 
a RE beam. The modeling of the current of this confined plasma will be discussed in 
Sec. [5] 

Monotonic radial profile of q(p ): The case g(0) > 1. In this case the m/n = 1/1 
mode does not play significant role in the formation of the stochastic zone in the plasma 
center. However, the m/n modes with n > 3 contribute much to the growth of stochastic 
zone and shrinkage of the intact magnetic surface p c . The examples of such stochastic 
magnetic fields are shown in Fig. [b] (a) and (b) corresponding to the values q((Y) = 1.1 
and q(0 ) = 1.2, respectively. 

The outward drift velocity v^r of such RE beams is significantly larger than the one 
in the case q( 0) < 1. This is because of the smaller RE beam radius p c and the higher 
toroidal electric field E v . Such RE beams decay in shorter times. 

Plasma with reversed magnetic shear. In this case the safety factor q(p) has 
a minimal value located at normalized radius p m /a ~ 0.4-0.6 and increases towards 
the center and the plasma edge. Figure [7] (a) and (b) show the example of the non¬ 
monotonic radial profile of the safety factor q{p) and the corresponding Poincare section 
of stochastic magnetic field lines (More detailed description of this case is given Sec. 7.3 of 
the Supplementary Part). The intact magnetic surface located near the shearless magnetic 
surface (red curve), i.e. the magnetic surface with a minimal value of the safety factor 
q(p), is not broken even at the relatively large magnetic perturbations. And it confines 
electrons in the central plasma region. 

Due to the relativ ely l arge confined area the RE beam would carry a large current 
Ip RE \ According to (3.2) the decay rate of this RE beam owing to the outward drift 
would be small. This effect probably explains the large RE current with a long lifetime 
observed in the TFTR tokamak during the disruption of plasmas with the reversed 
magnetic shear (Fredrickson et al|2015). 







Figure 7. (a) Radial profiles of the safety factor q{p ) in the plasma with the reversed magnetic 
shear, (b) Poincare sections of magnetic field lines in a pre-disruption plasma caused by several 
MHD modes. Red curve corresponds to the shearless magnetic surface. 


(RE) 


It follows from the 


3.3. Experimental evidences 

Existence of the finite interval of the RE currents Ij 

conjecture above that the RE current Ip RE is mainly determined by the current dis¬ 
tribution I p (p) in the predisruption plasma confined by the intact magnetic surface p c , 


i.e., Ip ~ Ip(p c ). Since p c is located between the magnetic surfaces p\ and P 3 corre¬ 
sponding to q(pi) = 1 and q(p 3 ) = 4/3, the RE current Ip RE> should be in a finite range. 
This expectation is in line with the experimental data presented in Fig. [4j One can see 
that the range of stable Ip RE ^ values shown in Fig. [d] b corresponds well to the space 
between resonant magnetic surfaces with q(p\) = 1 and q{p 3 ) = 4/3 [or g(p 2 ) = 3/2]. In 
Fig. [8]the radial profile of the pre-disruption plasma current I p (p) and the corresponding 

safety factor profile q(p) are plotted. Also values of Ip RE ' > found in other experiments on 
TEXTOR, see Ref. (Z eng et aL| |2013), lie in the same range. 

Since the g-value the plasma axis p = 0 is one of the major causes for uncertainties in 
the q (p)-profile, in Fig. [8]we show I and q profiles for g(0) = 0.75 and q( 0) = 0.8. These 
are in the range of g(0) values experimentally measured between sawtooth crashes in the 


TEXTOR tokamak (Soltwisch & Stodiek 1987 Soltwisch et al. 1987) (see also( Wesson 


2004), page 372). The values of g(0) measured after pellet injection in the DIII-D tokamak 
experiments are also close to these values (Izzo et al. 2012). Thus, small changes in q{ 0) 
still keeps the RE currents Ip RE ■* in the interval p\ < p < p 3 . The highest and lowest 
values of Ip RE ^ shown in Fig. [s] corresponding to the discharges #117859 and #120140, 
respectively, lie at the border of region pi < p < P 3 , p- 2 - They have the shortest duration 
time for the RE current decay (see Tableland Figs, [l] (b), [4j) . The presence of a several 
low-order m/n = 4/3, m/n = 3/2, and m/n = 1/1 resonant magnetic surfaces within 
the RE beam may lead to excitations of the corresponding MHD modes. The interactions 
of these modes may lead to the quick loss of REs due to the formation of stochastic zone 
at the edge of the RE beam (see Sec. [7|. 

Dependence on the level of magnetic perturbations. The existence of the intact 
magnetic surface p c between the q = 1 and q = 5/4 (or q = 4/3) rational magnetic 
surfaces and its location depends on the level magnetic perturbation cmhd (more exactly 
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Figure 8. Radial profile of the plasma current I p (p) ( |3.1| ) (solid curves 1 on l.h.s. axis) and the 
corresponding safety factor profile q(p) (dotted curves 2 on r.h.s. axis). The rectangular (red) 
dots correspond to the experimentally measured values of I i p RE ' 1 for several TEXTOR discharges. 
The plasma parameters are I p = 350 kA, Bo = 2.4 T, Ro = 1.75 m , a = 0.46 m. The values of 
qo = q( 0) are 0.75 and 0.8, respectively. The radii pi, p 2 , and po are the positions of the rational 
magnetic surfaces q(pi)=l, g(p2)=3/2, and g(p 3 )=4/3, respectively. 



on the spectrum B mn ). With increase of £mhd the radius p c shrinks and it can be 
broken at the certain critical perturbation level e cr . It leads to the total destruction of 
the confinement of plasma particles. This is in agreement with experimental observations 
on the existence of a critical magnetic perturbation level above which the runaway beams 
are not formed ( Zeng et q/.| 2013). 

The shrinkage of p c with increasing the magnetic perturbation cmhd leads to the 
decrease of the RE current I p RE ^ since Ip RE ^ ss I p (p c ). On the other hand if one as¬ 
sumes that the plasma current decay is caused by the radial transport of particles in the 
stochastic magnetic field, its decay rate dl p /dt should be proportional to t he square of 
the magnetic perturbation level cmhd, \{dl p /dt)\ oc | £mhd\~ (see Sec. 4.2). Therefore, 
one expects that to the higher values of | ( dl p / dt) | correspond the lower values of the RE 
current Ip RE \ This expectation is in agreement with the experimental values of these 
quantities presented in Fig. [4] (a). 

Synchrotron radiation pattern. The formation of the RE beam inside the intact 
magnetic surface can be also confirmed by the spatial profiles of the synchrotron radiation 
of high-energy REs with energies exceeding 25 MeV. Figure [9] shows the radial profiles 
of infrared radiation of the REs at the equatorial plane z = 0 for the two TEXTOR 
discharges. One can see that radiation is localized inside finite radial extent corresponding 
to the central region of plasma within the q = 4/3 magnetic surface (see Fig. [5](b)). One 
should note the radiation from the plasma edge regions, 15 < R < 1.6 and 2.1 < R < 2.2 
is due to thermal radiation of the wall elements. The outward shift of the radiation pattern 


is explained to some extend by the drift of RE beams discussed in Sec. 6.1 


Another indication of the formation of confined plasma beam is the rise of the tem¬ 
perature at the initial stage of the beam formation as seen in the ECE signals shown in 
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R [m] 


Figure 9. Radial profiles of the synchrotron radiation at the equatorial plane z = 0: curve 
1 corresponds to the discharge #117507 at the time instant t = 2.034324 s and curve 2 — to 
#120134 at t = 2.050284 s. 

Figs. |T] (a) and(b). It may occur due to the Ohmic heating of confined plasma by the 
induced toroidal electric field or by superthermal emission from high energy electrons. 
As thermal electrons are converted into runaway ones the beam temperature goes down. 


4. Thermal and current quench stages 

Our analysis in section 2.1 reveals that both the processes induced directly by injected 
atoms and strong radial transport along stochastic magnetic field lines created by MHD- 
perturbations can lead to heat losses from the plasma in the TQ stage of disruption. 
During this stage the temperature drops on a timescale of several hundreds of microsec¬ 
onds; current decay time is of (4 4- 6) ms in RE-free discharges and increases up to 0.1 s 
in discharges with RE generations. 

As one can see in Fig.l the CQ stage before the formation of RE beams coincides 
well with the time interval where magnetic perturbations are significant and the particle 
transport in stochastic magnetic field leads to the current decay. To study these processes 
we use the models for stochastic magnetic field and collisional transport of test particle 
described in Sec. 7 of Supplementary part. 

Figures 10 (a) and (b) shows the typical Poincare sections of field lines of this model 
in the runaway-free disruption case (a) and the case with RE generation (b). The per¬ 
turbation amplitudes e mn of all MHD modes, except (m = l,n = 1) mode, correspond 
to the twice larger value of cmhd than the case shown in Fig. [5j For the (m = 1, n = 1) 
mode e mn corresponds to the same value of cmhd- The relation between e mn and cmhd 
is Cmn — eMHDbmn/'&a, where is the toroidal magnetic flux at the plasma edge (see 
Sec. 7 of Supplementary part). 

In general the transport of heat and particles in the presence of R.MPs is a three- 
dimensional problem. Particularly, a stochastic magnetic field with the topological struc¬ 
tures like ones in Figs. [To] leads to poloidally and toroidally localized heat and particle 
deposition patterns on wall (Kruger et al. 2005). This is a general feature of open chaotic 


systems which has been observed in ergodic divertor tokamaks (see, e.g., (Finken et al. 


2005 Jakubowski et al. 2006 Abdullaev 2014)). The problem can be simplified when 


we are interested only in radial transport rate. It can be done by introducing the radial 
diffusion coefficient averaged over a poloidal angle. 
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Figure 10. Poincare sections of field lines in a pre-disruption plasma caused by several MHD 
modes: (a) runaway-free discharges; (b) with runaway electrons. The dimensionless MHD mode 
amplitudes are e mn = 8.68 x 10” 3 fe m „ with 6 mn = 1 for all modes in (a), and 6n = 1/4, b mn = 1 
for all other modes (n = 1, 2, m = 1 — 5) in (b). The safety factor at the magnetic axis is 
g(0) = 0.8 and at the plasma edge q a = 4.7. 


4.1. Heat transport 


The electron heat conductivity in a stochastic magnetic field has been assessed by di¬ 
verse approaches. We apply here the following formula for the electron heat diffusion \r 
deduced on the basis of simulations for transport of test particles, by taking into account 
coulomb collisions with background plasma species (Abdullaev 2013) (see also Sec. 10.4 
in ( Abdullaev||2014 )): 

Xr(p,T.) = (4.1) 

-L i Amfp 


where dm « Vp c = 1.33x 10 7 T^ 2 is the thermal velocity of electrons, Dpp(p) the diffusion 
coefficient of field lines (Dfl(p) ~ 10 -5 4- 10” 4 m), A m / p = 8.5 x 10 21 T 2 (p)/n(p) the 
mean free path length of electrons with the temperature T e and density n{p) measured 
in keV and m , respectively, and L c « 7rg(p)f?o is the characteristic connection length. 

A characteristic heat diffusion time one can estimate as Tp = a 1 /2\ r , where for \ r 
we assume its magnitude at the radial position p = 0.566 a. Before the disruption the 
local temperature here is of 0.6 keV. This provides \r = 287 m 2 /s and th = 3.68 x 10 -4 
s, i.e. of the order of the experimentally observed time for the plasma temperature drop 
during the TQ after disruption. 

For a quantitative analysis we have modeled the time evolution of the radial profile 
for the electron temperature averaged over the poloidal 9 and toroidal ip angles, T(p,t)- 
This is done by solving numerically the following diffusion equation: 


dT 

dt 


1 d 

pdp 


PXr(P,T) 


dT' 

dp 


(4.2) 


where the heat diffusivity is given by equation (4.1) and the applied boundary conditions 
are: dT(p)/dp = 0 at p = 0 and dT(p)/dp = —T/Sp at the plasma edge p = a, where 
Sp — 0.1 to is the characteristic e-folding length for the temperature decay in the scrape- 
off layer. 
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Figure 1 1. R adial profiles of the electron heat conductivity \ r {p,t) computed according to 
equation (4.11, (a), and of the electron temperature av erag ed over toroidal and poloidal angles, 
T e (p,t), found by solving heat conduction equation (4.21 numerically, (b), at different time 
moments after the disruption initiation. 


Below we consider an example of heat transport in a fully chaotic magnetic field shown 
in Fig. 10 (a). Figures 11 (a) and (b) show the radial profiles of the heat conductivity 
and the temperature at different times. One can see in this case the temperature drops 
and almost flattens within a time interval of order of 0.5 ms. 

In the situation with partially stochastic magnetic field, see Fig. 10 (b), anomalous 
turbulent transport in the very plasma core, p < 0.3, with intact magnetic surfaces is by 
two orders of magnitude smaller than in the outer region. In this case the temperature 
drop in the central plasma region could be explained by the effect of parallel transport 
in the magnetic island created by the m/n = 1/1 MHD mode and some level of stochas- 
tization in a thin layer near the island separatrix (see a review by Schuller (1995) for 
more details). 


4.2. Current quench stage 


As was discussed in Sec. 3.3 the current decay rate ( dl p /dt ) depends on the initial RE 
current I p RE ^ which in turn depends on the level of magnetic perturbations. Therefore, 
one can assume that the current decay rate ( dl p /dt) may implicitly depend on the level 
of magnetic perturbations. This assumption, however, may contradict to the traditional 
view that the current decay rate is determined by the time tqq = R/L, i.e., by the ratio 
of the plasma resistivity R to it’s inductance L: dlp/dt = —I p /tcq- In the simplest case 
when R and L are constants it leads to the exponential decay law I p (t) oc exp (—t/rcc))- 
However, in many tokamaks, particularly in the TEXTOR tokamak, the current decay 
evolution is better fitted by the linear function I p = a —bt rather than by the exponential 
function (see Fig. [3]). At present the reason of such a dependence is not quite clear. So 


other mechanisms may also play role in the CQ (see, e.g., Kadomtsev (1984); Gerhardt 


et al. (2009); Shibata et al. (2010) and references therein). 


Particularly, one cannot exclude that the poloidal and toroidal variations of the plasma 
current imposed by the initial MHD modes in the TQ stage do not disappear immediately 
with the temperature drop (see also Sec. [5] and Eqs. (5.1) and (5.51). This leads to the 
corresponding variations of the poloidal magnetic field which acts as non-axisymmetric 
magnetic perturbations. Below we discuss a possible role of the radial transport of par¬ 
ticles in a stochastic magnetic field created by such magnetic perturbations. 

The magnetic field structure before the CQ has been assumed similar to the one shown 
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Ti [keV] 

D p [m 2 /s] 

0.005 

0.0986057 

0.050 

0.386249 

0.100 

1.01251 

0.500 

6.46228 

1.000 

9.51915 

2.000 

13.1030 

4.000 

17.8366 

5.000 

23.7424 

10.00 

27.0265 


tcq = a 2 /2D P [s] 

1.072 

2.739 xlO' 1 
1.045 xlO' 1 
1.637 xlO~ 2 
1.111 xlO~ 2 
8.074 xl0~ 3 
5.932 xl0~ 3 
4.456 xl0~ 3 
3.915 xl0~ 3 


Table 2. Ambipolar diffusion coefficients D v of particles and the diffusion times tcq = a 2 /2D P 
from the stochastic zone at the different effective plasma temperatures. The plasma radius 
a = 0.46 m. 


in Fig. 10 The level of magnetic perturbations may be different from that during and 


after the thermal quench when strong magnetic fluctuations are present, see Fig. 0(d)- 
The timescale of the current decay is determined by the rate of radial particle transport 
in a stochastic magnetic field. This process has the ambipolar nature and it is strongly 
collisional due to the low plasma temperature. On other hand one expects that the 
toroidal electric field induced by the current decay also strongly affects the particle 
transport. Below we give a rough estimation of the particle transport rate based on 
the collisional test particle transport model. 


In Table [2] we have listed the ambipolar diffusion coefficients D p and the characteristic 
diffusion times tcq of particles at the different plasma temperatures in a stochastic 
magnetic field shown in Fig. [lOj The typical plasma temperature after the TQ is about 
from 5 eV to 50 eV. The average particle confinement time tcq at this temperature 
changes from 1 s to 0.3 s. These time scales are shortened if the magnetic perturbation 
level e 2 is larger. Since the diffusion coefficient D p oc e 2 and therefore tcq oc e -2 , then 
tcq can be reduced to one order smaller value for three times larger perturbation than 
in Fig. |10| This timescale is still much longer than the experimentally observed values. 
However, this collisional model does not takes account the effect of the inductive toroidal 
electric field. One expects that the acceleration of electrons and ions by the electric 
field increases the radial transport of particles. To include this effect in the collisional 
model one can assume that the effective temperature of the plasma is higher than the 
measured one. The particle diffusion time tcq for the effective temperature 2 keV is 
about 8x 10 -3 s. This timescale gives the average current decay rate \dl p /dt\ ss I p /tcq = 
0.35/(8.0 x 10 -3 ) ~ 44.0 MA/s which is of the order of the experimental measured rates 
given in Table [l] 


More rigorous approach to the particle transport in a stochastic magnetic field during 
the current decay stage would require a three-dimensional treatment of the problem. It 


should take into account not only the formation of the ambipolar electric potential (Spizzo 
et aZ. 2014) but also the inductive toroidal electric field which accompanies the process. 


The latter may lead to the directionality of the particle transport that eventually may 
influence on the random vertical displacement of runaway beams. Of course, the study 
of these complicated processes is out of the scope of the present work. 
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5. Modeling of post—disruption plasma 

The described scenario of plasma disruption with a RE beam allows one to model 
a post-disruption plasma. After establishing the runaway beam the current is localized 
inside the area enclosed by the last intact magnetic surface. In general the distribution 
of the current density j would depend not only on the radial coordinate p but also vary 
along the poloidal 6 and the toroidal ip angles due to the presence of the (m/n = 1/1) 
magnetic island. Such a post-disruption plasma current can be presented as a sum of 
two parts, 

j(p,Q,<P) = 3o(p)+ji(p,0,<p), (5-1) 

where jo(p) is the current density depending only on the radial coordinate p, and ji(p, 0 , tp) 
is the helical current which is a periodic function of the poloidal 6 and toroidal <p angles. 

The radial dependence of jo(p) can be also modeled by assuming that after the dis¬ 
ruption the current is uniformly distributed over confined area with the steep gradient at 
the beam edge p c . Calculations show that electron orbits does not significantly depend 
on the specific of the radial profile of jo(p). For our calculations of GC orbits we choose 
the following profile 

{ J 0 tanh \(p 2 c - p 2 )/A a ] , for p < p c , 

3o(p) = { ( 5 - 2 ) 

l 0, for p > p c , 


where Jo is the constant determined by the full current of the beam I p RE ^ and A a is 
the steepness parameter. The current flowing inside magnetic surface p, i.e., I p (p) = 

p 

2ttJo / jo(p l )p l dp l is given by 
o 


Ip(p) — ' 

I(H E) 

i In cosh (fp 2 — p 2 ]/ Aa ) 
In cosh [pi 2 / A a ) 


{ li RE \ 

for p > p c , 


(5.3) 


where Ip RE> is the full current of the confined area. 

One should also note the fact that after the TQ the plasma beam is shifted inwardly 
because of drop of plasma pressure. In the modeling this fact can be taken into account 
by assuming that the radial position of the center R a of the post-disruption plasma is 
different from the one of the pre-disruption plasma. The safety factor of the corresponding 
plasma is then given by 


l{p) = Qcyl(p)C{p/R a ), 


Qcyl(p) 


2np 2 Bo 

poRadpip) 


(5.4) 


where q cy i{p) is the safety factor of the cylindrical plasma, the function C(x) = 1 + 
A\X + A 2 X 2 + ■ ■ ■ is a function which takes into account the toroidicity of plasma. The 
coefficients A,, (i = 1,2 ,...) depends on the plasma pressure (Abdullaev et al. 1999 
Abdullaev||2006| . 


Figure 12 shows the radial profiles of I p (p) (solid curves 1-3 on the l.h.s. axis) (5.3) and 


the safety factor (5.41 (dashed curves 1' — 3' on the r.h.s. axis) for the three discharge 
parameter, respectively. Solid black curve 4 corresponds to the pre-disruption plasma 
current profile. We set the toroidal field magnitude Bq = 2.4 T, the beam center at 
R a = 1.7 m. The plasma radius a is found from the condition I p o(a) = Ip RE ' > where 
Ipoi r ) is the current profile of the pre-disruption plasma. The vertical dashed color 
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Figure 12. Radial profiles o f the plasma current I p (p) (5.31 (solid curves 1—3 oil l.h.s. axis), 


the safety factor profiles q{p) (5.41 (dashed curves 1' — 3' on r.h.s. axis), and curve 4 corresponds 

ARE) _ 


to the pre-disruption plasma current. The red curves 1 and 1’ correspond to I p 

ARE) 


165 kA, 

blue curves 2 and 2’ correspond to I p txI1 ‘ l = 230 kA, and green curves 3 and 3’ correspond to 
I p RE ' ) = 300 kA. The vertical solid arrows indicate the radii of plasma beam a, the vertical 
dashed arrows indicate the positions of resonant magnetic surfaces q = 1 and q = 3/2. The 


toroidal magnetic field B t = 2.4 T 
and the radius ao = 0.46 m. 


R a = 1.7 m, the pre-disruption plasma current I p = 350 kA 


arrows show the radial positions of the q = 1, q = 4/3, and q = 3/2 magnetic surfaces 
and the vertical solid arrows indicate the plasma radii a. 

Note that that the red curves 1 and 1' and green curves 3 and 3' in Fig. [Tg] correspond 
to the discharges with the lowest and highest values of Ip RE ^ shown in Fig.Tlbb). For the 
lowest value of Ip RE ' the radial position of the q = 1 magnetic surface is very close to 
the RE beam radius p c . For the highest value of Ip RE ' > the magnetic surfaces with q = 1, 
q = 4/3, and q = 3/2 are located inside the plasma region p < a. However, the radial 
position of the magnetic surface q = 3/2 is at the plasma edge. For the typical discharges 
like the one shown by blue curves the magnetic surface q = 1 is located relatively far 
from the plasma edge. 

The Fourier expansion of the helical current, ji(p, 9, ip), 

j 1 (p,9,ip) ='^2,j mn (p)cos(m8-rup + (t> mn ), (5.5) 

m,n 


is mainly dominated by the m/n — 1/1 component. This assumption is based on the 
analysis of numerous disruptions in the JET tokamak ( Gerasimov et ar§ 2014). 

We should assume that the value of the safety factor at the beam axis g(0) is less than 
unity. This assumption is supported by a number of experimental measurements of the 
current profile after the sawtooth crashes in the TEXTOR, the TFTR, and JET tokamaks 
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This model of the post-disruption plasma current describes only the initial stage of 
RE beam. During acceleration of electrons in the toroidal electric the RE orbits drift 
outward and their form evolves from the circular one to oval one. This process changes 
in turn the RE beam form and its current. The self-consistent description of the time 
evolution of the RE beam is a difficult problem. It is beyond of the scope of this present 
study. 


6. Evolution of GC orbits during acceleration 

Now we discuss the dynamics of RE orbits during the acceleration of electrons induced 
by the toroidal electric field in a toroidal post-disruption plasma. Mainly we study the 
outward drift of RE orbits and it’s role in the RE losses. However, we will not consider 
the processes of the generation and the proliferation of the RE population, the problems 
of stability of RE beams and related issues. These problems have been much discussed 
in literature. 


6.1. Outward drift of RE orbits 

First we consider the case of the axisymmetric plasma beam neglecting the helical mag¬ 
netic perturbations. The inductive toroidal electric field generated due to the current 
decay during the plasma disruption accelerates thermal electrons. This is an adiabatic 
process since the characteristic time of significant variation of energy is much larger than 
the transit time of electrons. Therefore the GC orbit slowly drifts outward without chang¬ 
ing the area of GC orbit in the poloidal plane which is an adiabatic invariant J or the 
action variable (see Sec. 6.1 of Supplementary part)). With increasing electron energy 
the topology of GC orbits also slowly changes from the circular one to the oval one. 
Starting from the certain critical energy E cr the adiabaticity of the process breaks and 
the GC orbit bifurcates by creating the unstable stagnation point (or X-points) inside 
the plasma region. With the further increase of energy the GC orbit crosses the separa- 
trix (a homoclinic orbit associated the X-point) and becomes unconfined. The value E cr 
depends on the plasma current I p . The described phenomenon is an addition mechanism 
of confinement loss of REs. Figure 13 (a) shows a typical evolution of a GC orbit in the 
presence of the toroidal electric field with the constant beam current I p = 100 kA and 
the loop voltage V = 40 V. 

One should note that the formation of the separatrix of RE GC orbits during the 
acceleration process in tokamaks has been first predicted in Ref. (Zehrfeld et al. 1981). 
The numerical study of this process in a realistic tokamak configuration has been carried 
out in Ref. (Wongrach et al. 20141. Particularly, it was shown that with increasing electron 
energy the area confined by the separatrix decreases and it vanishes when the energy 
exceeds a certain critical value £ cr , i.e. such electrons cannot be confined. The critical 
energy £ cr is proportional to the square root of the plasma current I p , £ cr oc y/lp- 

The described evolution of RE orbits is in agreement with the experimental observa¬ 
tion of the IR radiation patterns observed in the experiment in the TEXTOR tokamak 
( Wongrach et al.|2014 ) and in the DIII-D tokamak ( |Hollmann et aL||2013 ). The observa¬ 
tions clearly show the evolution of the spatial form of RE beam from crescent ones into 
oval ones with increasing the electron energies. 

The example of the time-evolution of GC orbits in the plasma beam with a time- 
varying current I p (t) and the loop voltage V(t) corresponding to the TEXTOR discharge 
#117527 is shown in Fig. 13 (b). To simplify the calculations of orbits we have assumed 
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Figure 13. (a) Evolution of the GC orbit of accelerating electron in the (R, Z )-plane at the 
constant plasma current I p = 100 kA in the presence of a constant toroidal electric field with 
V = 40 V. Dashed curve corresponds to the separatrix of GC orbits of electron of energy 
E = 27 MeV. (b) The same as in (a) but for the two GC orbits of accelerating electrons for the 
time- varying current I v (t) and the loop voltage V(t) corresponding to the TEXTOR discharge 
#117527. Green curves correspond to the separatrices of GC orbits. 


that the loop voltage V(t) is uniform in the poloidal section, i.e., it does not depend on 
the radial coordinate r and equal to the experimentally measured value at the limiter. 
However, this assumption only approximately describe the situation. To find more the 
exact magnitudes of the toroidal electric field during the runaway current decay one 
should solve the corresponding Maxwell equations. 

One of the important parameter of the GC orbit is the effective safety factor g e // 
defined as a ratio q e ff = Aip/2n where Ac p is the increment of the toroidal angle A per 
one poloidal turn. It is a function of the action variable J and particle energy E. For 
low-energy electrons the quantity q e ff(J,E) coincides with the safety factor q(p) of the 
equilibrium magnetic field. With increasing the electron energy the effective safety factor 
strongly deviates from q{p). With approaching RE energy E to the critical one E cr it 
diverges as 


q e ff{J , E) oc - In |E - E c , 


( 6 . 1 ) 


Figure 14 shows the typical time evolutions of the effective safety factors q e ff of two 
GC orbits during the electron acceleration in the conditions of the TEXTOR discharge 
#117527. 

Figure [is] shows the time-evolution of the outward drift velocity Vdr calculated nu¬ 
merically for the three different RE beam currents. It is quite well described by formula 
derived in Ref. (Abdullaev|[2015)) (see also Sec. 4 of Supplementary part) 


_ Ro-Ey (RT av \ 
Vdr ~ RB* V 1 ~ ~RqT ) ’ 


( 6 . 2 ) 


where B* = B z + F(E) is the effective poloidal magnetic field, B z is the 2 — component 
of the poloidal magnetic field at the equatorial plane 2 = 0, F(E) is the term depending 
on a particle energy. The quantity 


T — 

av — 


27T q e f jRq 


(6.3) 


is the average transit time, E v is the toroidal electric field strength, T is the transit time 
of orbit, v v is the toroidal velocity. 

The expression (6.2) describes the creation of the X-point and the separatrix of RE 
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Figure 14. Time-evolution of the effective safety factors q e ff (l.h.s. axis) and electron energies 
(r.h.s. axis) during the acceleration in the discharge #117527. Blue curves correspond to the orbit 
launched at the coordinate (R = 160, Z = 0) cm, and red curves to the one with (R = 165, Z = 0) 
cm. Horizontal lines correspond to q(t) = m where m = 1, 2,... are the integer numbers. 



Figure 15. Drift velocities of innermost Vdf(Ri) (curves 1, 2, and 3) and outermost Vdf(Ro) 
(curves T, 2', and 3') points of orbits for the different plasma current: curves 1 and 1' correspond 
to the plasma current I p = 100 kA, curves 2 and 2' correspond to I p = 200 kA, and curves 3 
and 3' correspond to I p = 300 kA. Curve 4 describes the increase of energy E (right hand axis). 
The toroidal field Bt = 2.5 T, major radius Ro = 175 cm, minor radius a =46 cm, the loop 
voltage Vioop = 5 V. Note that Vdr is multiplied to the proportionality factor /. 


orbits at the critical energy E cr . This phenomenon is related with the appearance of 
zeroes of the effective poloidal magnetic field B* z at E = E cr and the certain radial 


distance R = R s within the plasma region (see Abdullaev |2015) for details). 
At \B Z \ \F(E)\, and T av « T the formula (6.21 is reduced to 


Vdr — 


qRf 

Bn 


(R — Ro)E lf 
RB Z 


(6.4) 


obtained by |Guan et al. (2010); Qin et al. (2011) for the circular orbits. Here q = 
(R — Rq)Bo/B z R is the safety factor of magnetic field. As seen from Fig. 15 the formulas 


(6.2) and (6.4) give the correct dependence of Vdr on the plasma current I p , Vdr oc / 


-l 


because B z oc I p . 


However, the formula (6.4) does not describe the situation when the GC orbits take an 


oval form with increasing the energy similar to the ones shown in Figs. 13 (a) and (b). 


From the latter it follows that the average outward velocity Vdr of the innermost part of 
the orbit is approximately equal 0.6 m/s and 8 m/s of the outermost part of the orbit. 
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6.2. RE current decay 

The rate dl p /dt of the runaway current loss due to described outward drift of orbits can 
be roughly estimated as follow. This loss mechanism is mainly caused by the shrinkage 
of the beam radius p c . The rate of such a shrinkage dp c /dt is of order of the average 
outward velocity Vd r ■ Since I p oc p 2 cl we have 


dip 

dt 


dl p dp c 
dp c dt 


2 I p E v 

- Vdr OC -. 

Pc Pc 


(6.5) 


For the typical values of I p ~ 0.2 MA, p c ~ 0.2 m, and v,j r ~ 1 m/s one has dl p /dt ~ 4 
MA/s. This estimation is of order of the experimentally measured average decay rate of 
the runaway current listed in Table [T] 

Since the safety factor q , as well as q e f / of RE beams is about unity, q e f / ~ 1 then 
the outward drift may slowdown for the higher values of the toroidal magnetic field Bq . 
However, much the higher toroidal electric field may compensate this effect so that 
the decay time of RE currents in large tokamaks, like ITER, may have the same order 
as in smaller tokamaks. 

One should also note that the outward drift velocity Vd r is proportional to the inverse 
aspect ratio of tokamaks, v,i r oc a/Ro (Abdullaev 2015). It means that the RE current 
loss due to the outward drift of orbits in spherical tokamaks would be larger than in 
standard tokamaks, so that the RE electron would cease faster. This effect could be one 


of reasons of the absence of REs during disruptions in NSTX tokamaks (Gerhardt et al. 


2009). 


Beside of outward orbit drifts the RE current losses is also caused by the internal MHD 
mode which be discussed in the next section. The collisions of REs with neutral particles 
may also contribute to the RE losses. 


7. Effect of magnetic perturbations 

The effect of the magnetic perturbations on electrons in the post-disruption current 
beam strongly depends on its safety factor profile q(p), the spectrum of magnetic per¬ 
turbations, and the electron energy. To explain this effect we consider the simplified 
version of GC motion equations in the presence of magnetic perturbations. (The rigorous 
consideration of this problem is given in Sec. 6 of Supplementary Part). 

The particle drift motion in the presence perturbations can be presented by Hamilto¬ 
nian equations similar to the equations for magnetic field lines, 

dd dK dJ dK 

dtp dJ ’ dp dd z ’ 

with the Hamiltonian K = K(d, J, E, p) with the canonical variables (d,J), and the 
toroidal p as the time-like variable. In the absence of perturbations GC orbits wound the 
drift surfaces J =const and the poloidal angle d is a linear function p, d = p/q e ff(J, E) + 
d. In the presence of perturbations Hamiltonian H can be presented as a sum 

K = j E) + eh\(d z , J, E, p). (7.2) 

Since the perturbation are periodic in poloidal and toroidal angles and in time it can be 
presented by a Fourier series 

Ki{d z , J, E,p) = Kmn{J , E) exp [i ( md - np)\. 


mn 


(7.3) 
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Figure 16. (a) Spectrum of perturbations K mn and (b) corresponding RE orbits with different 
energies E. Curves 1-7 correspond to RE energies 10 keV, 20 MeV, 30 MeV, 40 MeV, 42 MeV, 
42.5 MeV, and 42.7 MeV, respectively. Curve 8 corresponds to the separatrix with the critical 
energy E cr = 42.646 MeV. The plasma current I p = 150 kA, the toroidal field Bo = 2.5 T. The 
toroidal mode number n = 1. 


The strongest influence of perturbation on particles takes place on the (to, n ) resonant 
drift surfaces, i.e., 


m = nq eff {J,E), 

originating from the (to, n ) term in ( |7.3| ) with the amplitude K„ 
termined by the magnetic perturbation spectrum 6 mn , 

r 2n 


(7.4) 

,.(,/, E). They are de- 


K r , 




mn I did exp [i (mid - to'i?m)] , 
o 


(7.5) 


where iJ m the poloidal angle associated with magnetic field lines is a function of id as 
well as particle energy E. 

For low-energy electrons (up to 5 MeV) the spectrum of amplitudes K mn (J , E ) weakly 
depends on energy E and close to the spectrum of magnetic perturbations b mn of 
(to, n)—th modes. With increasing the energy the spectrum of perturbations K mn (J,E ) 
deviates from b mn and acquires more higher poloidal harmonics to. The example of the 
poloidal spectra of perturbation K mn (J,E) for different particle energies is shown in 


Fig. 16 (a). The corresponding unperturbed orbits are plotted in Fig. [16] (b). It is as¬ 
sumed that the magnetic perturbation contains a single (to = l,n = 1) mode. For the 
low energy electrons with E < 10 MeV the spectrum K rnn contains the predominant 
to = 1 mode. 

With increasing the energy the amplitudes K m n of higher to also grow and the width 
of the poloidal spectrum /\ mra in to becomes wider as shown in Fig. 16 (b). For the 
spectrum K mn one can obtain the following asymptotical formula for the orbits close to 
the separatrix (see Sec. 3.4 in (Abdullaev|2014)) 

1 ( rnC\ 

exp ( — -—), (7.6) 


K„ 


deff 


leff 


where C is a finite constant, and the effective safety factor q e ff diverges as ( |6.1[ ). 

As was shown in Sec. [5] (see also Fig. 12) the typical values of q(p) varies between 
q( 0) ~ 0.7 4- 0.8 at the magnetic axis and q(a) < 1.5 at the plasma edge. Therefore, the 
strongest effect of the RMPs on electron orbits may expect if its spectrum b mn contains a 
sufficient number of (to, n)— components that are resonant to the magnetic surfaces with 
q in the interval g(0) < q = m/n < q(a) that would create a stochastic zone of magnetic 
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Figure 17. Poincare sections of RE orbits in the ( R , Z)-plane: Red dots correspond to RE orbits 
with energies E = 1 MeV (a) and E = 20 MeV (b) , blue dots correspond magnetic field lines 
(a); green curve is the separatrix. The perturbation parameter e = 10 —5 , the plasma currents 
I p = 200 kA (a) and I p = 100 kA (b). 




field lines. The electrons from this stochastic layer would be then radially transported to 
wall. 

Below we discuss the influence of magnetic perturbations on RE orbits for the two 
specific cases. First we consider the effect of internal single helical magnetic field, and 
then we analyze the effect of the external RMPs, namely, the TEXTOR-DED on the 
confinement of REs. 


7.1. Effect of a single helical magnetic field 


Assume the magnetic perturbation (3.3) contains the single (m = l,n = 1) MHD mode 
as was proposed in the model of the post-disruption current beam described in Sec. [5j 
For the low-energy electrons it creates a single island structure since the deviations of 
their GC orbits from the magnetic surfaces is small. Such a system is stable because the 
single MHD mode does not create stochasticity of magnetic field lines. The example of 
this case is shown in Fig. [l7] (a) by the Poincare sections of RE orbits (red dots) and 
magnetic field lines (blue dots). 

With increasing the energy of electrons and decreasing the beam current the electron’s 
GC orbits strongly deviate from the magnetic field lines. The effective safety factor q e ff 
of the GC orbit increases as the electron energy grows as was shown in Fig.[l4j At certain 
time instants the value of g e // reaches the integer value so that the resonant condition 
may be satisfies for the higher harmonics (m > 1, n > 1) of the GC orbits with the 
(m = 1 ,n= 1) magnetic perturbation. This generates a number of island chains of GC 
orbits. The interaction of several such island structures may even lead to the formation 
of the stochastic layer near the separatrices [see Figs. 13 (a) and (b)]. 

Figure 17 (b) illustrates the typical structure of high-energy electrons in the presence 
of the internal helical magnetic field with a single (m = 1, n = 1) mode. Such a structure 
leads to the widening area of lost electrons and decreasing the critical energy E cr . The 
characteristic escape time of REs from the stochastic layer is of order of 10 /us. Sudden 
RE bursts in many discharges is probably related with the loss REs from the stochastic 
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layer. Occurrence of the MHD mode signals accompanied these events will be discussed 
in the next Sec. 17.31 

As was discussed above in Secs. [2] and [5] (also Figs. [T] (b) and 12) there are some 
exceptional discharges (for example #117859) with the highest RE current and several 
low-order rational surfaces within plasma beam. Such a beam can be easily destabilized 
by the magnetic perturbations containing several MHD modes with low-order (m, n) 
numbers. Such a magnetic perturbation may affect strongly on electrons creating the 
chaotic zone at the beam edge open to wall. Such an effect probably explains the sudden 
lost of REs at certain times seen in Fig. 0 GO- 


7.2. Influence of the TEXTOR-DED 

The coil configuration of the TEXTOR DED is designed to have the poloidal spectra 
of magnetic perturbations localized near the magnetic surface q = 3 of the flat-top 
plasma discharges (see Sec. 5.2 of Supplementary part). Therefore, these perturbations 
do not contain a necessary number of resonant components to create a stochastic zone 
of magnetic field lines in the post-disruption current beam with the safety factor q lying 
between g(0) < 1 and q(a) < 1.5. 

In the so-called 3/1 operational mode with the predominant toroidal mode n = 1 
toroidal there is only one (m = 1 ,n = 1) component resonant to the magnetic surface 
q = 1. The similar situation takes place in the 6/2 mode (n = 2) with the resonant 
component (m = 2, n = 2). [There are no magnetic surfaces in the plasma region that 
are resonant to the components (m = l,n = 2) and (m = 3 ,n = 2)]. On the other 
hand this resonant component of the DED field is weak since it is located away from 
the maximum of the spectrum. Therefore, the effect of the DED on the RE beam does 
not create the stochastic zone of magnetic field lines from which electrons would escape 
to wall as in the case of the stochastic zone in a flat-top plasma operation. The only 
m/n = 1/1 component of the DED perturbations may create an island structure near 
the q = 1 magnetic surface similar to one shown in Fig. |Tt] (a). 

With increasing the energy of REs and decreasing the plasma current the DED per¬ 
turbation starts to affect on REs because of appearance of high-mode resonances q e f / = 
m/n similar to the case discussed in Sec. 7.1 It generates the structures with islands 
and a stochastic layer. Figures 18 (a) and (b) show the typical Poincare sections of GC 
orbits of energetic electrons affected by the TEXTOR-DED: (a) corresponds to the 3/1 
mode with the DED current Id e d = 3 kA; (b) corresponds to the 6/2 mode with I ded = 7 
kA. The particle energy is taken E = 20 MeV, the plasma current I p = 94 kA, and the 
toroidal field B$ = 2.5 T. These structures explain the fast decay of RE current in its 
final stage accompanied by spikes in the scintillation probe (see Fig. 0(b)). 

The structures of RE orbits shown in Figs. 17 (b) and 18 correspond to the final 
termination stages of RE current. They have features which are characteristic for the 
so-called stable and unstable manifolds created by the splitting of separatrices (see, 
e.g., (Abdullaev 2014)). They lead to the toroidally and poloidally localized deposition 


patterns of REs on wall. Toroidal peaking and spatial-temporal evolution of hard X-ray 


emission in the final stage of RE current loss observed in DIII-D experiments (James 
et al. |[2012) is consistent the described topology of REs. 


The experimental observations in the TEXTOR-DED have indeed showed that the 
RMPs field which switched on just after the TQ does not affect on the radial transport 


and the loss of low-energy electrons (Koslowski et al. 2014 Wongrach et al. 2015) (see 
also Table [I]). This is mainly because of the mentioned features of the poloidal and 
toroidal spectra of the DED field. 
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Figure 18. Poincare sections of RE orbits in the ( R , Z )-plane of the RE orbits of energy E = 20 
MeV in the presence of the magnetic perturbations of the TEXTOR-DED. (a) corresponds to 
the n = 1 mode and the DED current Ided = 3 kA; (b) corresponds to the n = 2 mode and 
Ided = 7 kA. The plasma current I p = 94 kA, the toroidal field Bq = 2.4 T. 


7.3. Generation of magnetic perturbations by high-energy electrons 
The above mentioned in Sec. [2] the occurrence of the MHD activities during the sudden 
RE bursts can be explained by the nonlinear interaction of high-energy electrons with 
the (to = l,n = 1) MHD mode. The MHD magnetic perturbations with mode numbers 
(to, n)— higher than the initial (to = l,n = 1) mode can be generated during the accel¬ 
eration process of the REs. At certain energy of REs their orbits strongly deviate from 
the magnetic surfaces which creates in turn higher (to, n)— harmonics, (to > 1 ,n > 1), 
of the MHD (to = l,n = 1) mode (7.5). The resonant interaction of RE orbits with 
these harmonics leads to the redistribution of corresponding current near these orbits 
according to the helicity of these modes. Therefore, the current density (5.51 acquires 
higher (to, n)— components j mn which in turn generates the corresponding MHD modes. 
The bursts of magnetic activities accompanied by accompanied by sudden runaway cur¬ 
rent drops observed in experiments (see Figs, [l] (a) and (b)) are probably related to the 
described phenomenon. 


8. Summary 

Based on the analysis of numerous experimental data obtained in the TEXTOR 
tokamak we have proposed a possible mechanism of the plasma disruption with the 
formation of RE beams. The plasma disruption starts due to a large-scale magnetic 
stochasticity caused by nonlinearly excited of MHD modes with low (to, n) numbers 
(■ m/n = 1/1, 2/1, 3/2, 5/2,...). The RE beam is formed in the central plasma region 
confined by the intact magnetic surface. Its location depends on the safety factor profile 
q(p) and the spectrum of MHD modes. In the cases of plasmas with the monotonic profile 
of q(p) and at the sufficiently small amplitude of the m/n = 1/1 mode the most stable 
RE beams are formed by the intact magnetic surface located between the magnetic sur¬ 
face (? = 1 and the closest low-order rational surface q = m/n > 1. Depending on the 
spectrum of magnetic perturbations this rational magnetic surface could be one of these 
ones: q = 4/3, q = 5/4 or q = 3/2. 

Such an outermost intact magnetic surface forms the transport barrier for particles in 
the central plasma region. Electrons in this confined region are accelerated by the indue- 
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tive toroidal electric field. Such a situation occurs, for instance, in plasma disruptions 
with runaway beams initiated by the Argon gas injection. Heavy Ar atoms do not pen¬ 
etrate sufficiently deep into the plasma and therefore they do not excite the m/n = 1/1 
mode with the amplitude necessary to create the fully chaotic magnetic field. On the 
other hand the injection of the lighter noble gases neon and helium does not generate 
runaways since the light gases penetrate deeper into the plasma and excite the large- 
amplitude ( m/n = 1/1) mode. 

During disruptions of tokamak plasmas with the reversed magnetic shear RE beams 
can be formed in the central plasma region confined by the shearless magnetic surface. 
The latter cannot be broken even at the relatively large magnetic perturbations and 
it acts as a robust transport barrier to the parallel motion of particles along chaotic 
magnetic field lines. One expects that electrons confined by this intact magnetic surface 
form a relatively stable RE beam with a large transversal size. Experimental observations 
of RE beams with long confinement times during disruptions of plasma with the reversed 
magnetic shear in the TFTR tokamak (Fredrickson et al. j 20151, probably, supports this 
expectation. 

Based on this scenario we proposed the models of the pre-disruption and post-disruption 
plasmas with REs to study the processes of thermal and currents quenches, the runaway 
current losses. The model of magnetic field was proposed to describe the large-scale 
magnetic stochasticity due to interaction of low-mode-number MHD modes. The radial 
transport of heat and particles in a stochastic magnetic field are studied using the colli- 
sional diffusional models. It was shown that the temperature drop during the fast phase 
of disruption is caused by the radial heat transport determined by the collisional electron 
transport in a stochastic magnetic field. We have estimated a current decay time using 
the ambipolar collisional particle transport model. The dynamics of RE orbits in a post¬ 
disruption plasma in the presence of the inductive toroidal electric field is investigated by 
integrating the equations of guiding center motion. We analyzed the effect of the internal 
MHD mode and external RMPs on the topology of RE orbits. 

The new model reproduces for the first time remarkably well the essential features of 
the measurements: 

(a) The outer part of the plasma is clearly ergodized while the inner section is still 
intact. This agrees with the observation that the runaways are only seen in the inner half 
of the torus while they are obviously quickly lost from the outer part. 

(b) one observes a short tiny spike during the energy quench; we have described this 
spike previously; this spike is attributed to the loss of runaways born at the start up of 
the discharge from the ergodic zone. 

(c) In the case of disruptions caused by injection of light He and Ne gases the energy 
quench is due to processes directly related to the penetration of neutral impurity into 
the plasma; if heavier argon is injected fast electron transport in stochastic magnetic 
field is of more importance. Stochastic motion of plasma particles is also responsible 
for the current decay; however due to the ambipolarity of particle losses the ion motion 
essentially affects the duration of the CQ stage. The estimations of the energy quench 
and the current decay times based on the models presented agree well with observations. 

(d) The slow decay of the RE current in the plateau phase is explained by the loss of 
runaways due to two effects: (*) an outward shift of the runaways due to their continuous 
acceleration and the subsequent loss at the wall; (ii) by the formation of a stochastic 
layer of high-energy REs at the beam edge in the presence of the m/n = 1/1 MHD 
mode. 

(e) The effect of the external resonant magnetic perturbations on low-energy electrons 
(up to 5-10 MeV) is weak and does not cause their loss. 
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The new mechanism explains well the observed disruptions in present day tokamaks. One 
can expect the following consequences, e.g., for ITER. 

(1) The structure of the stochastic zone during the TQ allows persistence of preexisting- 
runaways through this phase such that they act as seeds during the following phase of 
high loop voltage. 

(2) The decay phase of the REs is rather long such that REs can acquire very high 
energy. 

(3) External magnetic perturbations acting on REs seems are little promising unless 
the core of a RE beam can be ergodized. 

(4) A means for eliminating the REs completely is the injection of about 10 25 molecules 
H 2 or D 2 into the discharge (Hender et al. 2007). This massive gas injection may impose 
a heavy load on the cryo-pumping system. 
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